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ABSTRACT 

We exploit a set of high signal-to-noise (~ 70), low-resolution (R ~ 800) quasar spectra to search for the signature of the so-called 
proximity effect in the Hi Lyo- forest. Our sample consists of 17 bright quasars in the redshift range 2.7 < z < 4.1. Analysing 
the spectra with the flux transmission technique, we detect the proximity effect in the sample at high significance. We use this to 
estimate the average intensity of the metagalactic UV background, assuming it to be constant over this redshift range. We obtain a 
value of y = (9 ± 4) X 10"^^ erg cm"^ s"' Hz"' sr"', in good agreement with previous measurements at similar z- We then apply the 
same procedure to individual lines of sight, finding a significant deficit in the effective optical depth close to the emission redshift 
in every single object except one (which by a different line of evidence does nevertheless show a noticeable proximity effect). Thus, 
we clearly see the proximity effect as a universal phenomenon associated with individual quasars. Using extensive Monte-Carlo 
simulations to quantify the error budget, we assess the expected statistical scatter in the strength of the proximity effect due to shot 
noise (cosmic variance). The observed scatter is larger than the predicted one, so that additional sources of scatter are required. We 
rule out a dispersion of spectral slopes as a significant contributor Possible effects are long time-scale variability of the quasars 
and/or gravitational clustering of Lyo- forest lines. We speculate on the possibility of using the proximity effect as a tool to constrain 
individual quasar ages, finding that ages between ~ 10* and ~ 10^ yrs might produce a characteristic signature in the optical depth 
profile towards the QSO. We identify one possible candidate for this effect in our sample. 

Key words, diffuse radiation ~ intergalactic medium ~ quasars: absorption lines 



1. Introduction 

The multitude of absorption lines seen in the spectra of high 
redshift quasars gives important information about the state of 
matter in the universe, tracing the physical conditions of the in- 
[ tergalactic medium (IGM) at various epochs. It is commonly be- 
lieved that for column densities up to A^hi ~ 10'^^ cm"^, the 
absorbers are in photoionisation equilibrium with a metagalac- 
tic ultraviolet background field (UVB), composed of the inte- 
gral over all sources of UV radiation (essentially, star-forming 
galaxies and quasars). High-resolution spectra of the Lyman for- 
est provide not only a rather detailed statistical characterisation 
of the absorber properties such as line number densities as well 
as temperature and density distributions (e.g.. lKim et al.ll2(j0l]) . 
but also phy sical parameters such as Hi and Hei i photoioni- 
sation rates jRauch et alj Il997t iFardal et all 1 19981) . which di- 
rectly relate to the intensity of the UVB, as described recently 
by Bolton et al. (2005) invoking hydrodynamical simulations. 
Independently, the UVB has been successfully synthesised by 
combining the observed quasar luminosity function and the UV 
emission from galaxies (although the latter is still very uncer- 
tain) with the propagation of diffuse radiation jHaardt & Madaul 
[1996.) . 

In the vicinity of strong UV sources such as bright quasars, 
the H I photoionisation rate should locally increase, further re- 
ducing the density of residual neutral hydrogen. This enhances 
the transparency of the IGM to Hi ionising radiation and 
should become observable as a weakening of the Lyman for- 
est absorption near such sources. Such an effect has first been 
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noted by Carswell et al.l (Il982l) and was lat er baptised 'Inverse' 
(lMurdochetal ."1986') or 'Proximity Effect' B aitlik et al.1 ( Il988l 
hereafter BD088). Its prime application has been so far the pos- 
sibility to derive an independent estimate of the UVB inten- 
sity, by measuring the reduction of column densities (against 
the global evolution of absorption line density increasing with 
redshift) and combining this with the QSO luminosity at the 
Lyman limit, assumed to be known. The best constraints on the 
UVB using the proximity e ffect stem from the combined analy- 
sis of large quasar samples (ICooke et al.lll997l;IScott et al.ll2000l; 
Liske & Williger 2001), yielding mostly values consistent with 
the above quoted other methods. However, the uncertainties are 
still substantial. Besides the problem that a limited number of 
lines of sight always suffers from 'cosmic variance', there may 
also be systematic biases. In particular, if QSOs reside in intrin- 
sically overdense environments then the sign ature of the proxim- 
ity ef f ect will appear weak er than predicted jLoeb & EisensteinI 
ll995URollindeetal.ll2005h . Another uncertainty is the possibly 
limited lifetime of quasars. On the other hand, the proximity ef- 
fect may also be used to derive constraints on this import ant, but 
largel y unknown astrophysical quantity (e.g., Pentericci et al.l 
l2002h . 

In this paper we present an exploitation of new observa- 
tional material in terms of the proximity effect (Sect. |2]i. Rather 
than the traditional line counting we use the more sensitive 
flux transmission statistic to search for proximity effect signa- 
tures, augmented by extensive Monte-Carlo simulations to cal- 
ibrate the systematic and statistical eiTors (Sect. O. We deliver 
our results in Sections 14.11 and |5] Firstly we briefly present an 
analysis of the combined sample of 17 QSO spectra and de- 
rive an estimate of the UV background intensity. Secondly we 
demonstrate that the effect is measurable on single sightlines 



2 



A. Dall'Aglio et al.: The line-of-sight proximity effect in individuai quasar spectra 



2 
1.5 

1 r 
0.5 


1.2 

1 

0.8 
0.6 1^ 
0.4 
0.2 







4500 
Wavelength [A] 



Fig. 1. Example of a quasar spectrum taken 
from the dataset (HE 0940-1050, z = 3.088). 
In the upper panel the flux density in units of 
erg cm"- s"' A"' as function of wavelength 
is presented with the uncorrected (dotted line), 
corrected (long-dashed line) continuum level 
and the noise of the spectrum (dashed-dotted 
line). In the lower panel the fitted transmission 
is shown in gray while the corrected continuum 
is plotted in black (See Sect. 12.1 1 and 13.2.31 for 
details). In both panels the vertical dotted line 
represent the emission redshift. 



(IWilliger et al.1 1 1 9941: iLu et Zl 1 1 9961: LSavaglio et all T997h. not 
only statistically in large samp les ( Baitlik et aljl988l: Scott et all 
l200(]tlLiske & Wimgerll200ll) . 

The effect can actually be systematically detected in each 
single QSO of our sample (with one special case that is discussed 
separately). We speculate about the possibilities to detect signa- 
tures of finite quasar ages by virtue of the proximity effect. 
Throughout this paper, we assume a flat A-Universe with 
= 0.3 and Q.a = 0.7 and Ho = Tlkms' Mpc 



0.1. For those spectra taken in clear to photometric condi- 
tions (6 in total), the correc t ed m agnitudes matched the values 
from IVeron-Cettv & Veroiil (1200 6) within SVmag < 0.1. The 
V magnitudes of the remaining objects, after slit loss correc- 
tion, w ere systematically lower than the Veron-Cettv & VeronI 
(|2006|) values by not more than 0.3 mag. Since those data 
were taken in relatively poo r sky conditions, we adopted the 
IVeron-Cettv & VeronI ( l2006l) values and associated larger uncer- 
tainties (^Vmag < 0.3). 



2. Data 

2.1. Observations and data reduction 

Our data were obtained at the ESO-VLT UT4 (Yepun) in service 
mode between Oct 2002 and Apr 2003. We used FORS2 in long 
slit spectroscopy mode with the 600B grism, covering the range 
33 15-6360 A. A longitudinal atmospheric dispersion corrector 
(LADC) was used to account for differential refraction effects. 
With a slit width of 1", the resolution is ~ 800 (4 A FWHM). 
Table [T] summarises the observations for the full sample. In to- 
tal we observed 17 QSOs, always fully covering the Ly Forest 
spectrum between the Lya and Ly/3 emission lines. 

The spectra were reduced using IRAF standard procedures. 
Each exposure was bias-corrected and flat-fielded, cosmic rays 
were marked using a k-ct clipping algorithm, and the images 
were finally sky-subtracted. The extracted spectra were cali- 
brated in wavelength and flux and corrected for vacuum and he- 
liocentric shifts. The galactic extinction was taken into account 
assuming the E( B - V) estima t ions b y iSchlegel et al.l (1 19981) . 
together with the Cardelli et al.l (1 19891) extinction curve assum- 
ing Ry = 3.1. The two exposures available for each target were 
coadded with inverse variance weights, yielding a typical signal- 
to-noise ratio (S/N) of ~ 70 in the Lyman forest. 

2.2. Quasar magnitudes 

An evaluation of the proximity effect relies on the accurate 
knowledge of quasar fluxes. Even though absolute spectropho- 
tometry is compromised by intrinsic quasar variability, we 
reached reasonable accuracy in almost all spectra (See Tab. [T] 
for details). We accounted for slit losses modelling a gaussian 
point spread function with FWHM given by the average see- 
ing during the observations. Centering the 1 arcsec slit on the 
Gaussian centroid, we corrected for the flux falling outside the 
slit in all our spectra. In addition to our data, we used images 
obtained at ESO-VLT UTl (Antu) in Nov 2004 under photo- 
metric conditions, covering 3 of our fields and SDSS infor- 
mation for 4 quasars yielding consistency to within ^Vmag < 



2.3. Systemic quasar redsliifts 

Our spectra cover a sufficient range in wavelength so that we 
could measure the redshift of each quasar from more than one 
emission line. All measured redshifts are compiled in Table|2] In 
order to adopt a sy stemic redshift we used low-ionis ation lines 
whenever possible (lGaskel]||l98"2t iTvtler & Fan|[T99 2). For ob- 
jects where this was unfeasible, we used the redshift from high- 
ionisation lines with a statistical correction, determined from 
the average shift between Siin-Oi and Siiv-i-Oiv]. The agree- 
ment between Si in-O i and C ii estimates is generally good, even 
though the second line is usually very weak and rather broad. 
For Q 0000-26 only the asymmetric Lya line with strong as- 
sociated absorption was covered b y our spectrum; for th is ob- 
ject we adopted the redshift from ISchneider et"an (I199L) . For 
Q 0347-383 o ur corrected va lue is in agreement with the esti- 
mation done bv lSteidell (119901) . 

2.4. Quasar continuum 

The analysis of absorption lines requires a normalisation to the 
QSO continuum. We explored two types of continuum estimates: 
(i) a global power law (/(v) oc v"), excluding emission and ab- 
sorption regions, used to estimate the quasar flux at the Lyman 
Umit; (ii) a more local estimate that also includes the broad emis- 
sion lines as quasi-continuum. For this tas k we developed an 
aut omatic algorithm, follo wing the work by I Young et al.l ( Il979h 
and ICarswell et al.l (Il982i) . which perform a cubic spline inter- 
polation based on adaptive intervals along the spectrum with re- 
spect to the continuum slope. The points for the spline interpo- 
lation were chosen starting from a regular sampling of the spec- 
trum with a binning that becomes finer whenever the slope of the 
computed continuum exceeds a given threshold. This is done in 
order to better reproduce the wings of emission lines. 

In Sect. 13.2.31 below we assess the expected errors (aris- 
ing mainly from line crowding) associated with this process. 
Figure [T] shows a sample quasar spectrum together with the esti- 
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QSO 


V mag 


Z m 


Exp. Time Seeing 
(s) (arcsec) 


Airmass 


Sky Condition^ 


Obs. Date 


Ref. Mag 


CTQ 0247 


17.4 


3.025 


750 


1.29 


1.07 


CL, WI 


Dec. 8, 2002 


4 


CTQ 1005 


18.4 


3.205 


1500 


1.13 


1.21 


TN 


Jan. 9, 2003 


3 


CTQ 0460 


17.5 


3.139 


900 


1.45 


1.10 


PH 


Dec. 23, 2002 


4 


H 0055-2659 


17.5 


3.665 


600 


1.34 


1.02 


CL, WI 


Dec. 8, 2002 


4 


HE 0940-1050 


16.4 


3.086 


600 


0.86 


1.12 


TN, TK 


Nov. 26, 2002 


3 


HE 2243-6031 


16.4 


3.010 


600 


1.02 


1.24 


CL, PH 


Nov. 9, 2002 


4 


HE 2347-4342 


16.7 


2.885 


600 


0.77 


1.06 


PH 


Nov. 10, 2002 


1 


PKS 2126-15 


17.0 


3.285 


600 


0.97 


1.20 


PH 


Oct. 29, 2002 


4 


Q 0000-26 


18.0 


4.098 


600 


1.31 


1.03 


TN, CL 


Nov. 7, 2002 


4 


n 0009 All 


1 7 1 


Z. / D / 


600 


1.39 


1 J 


TIC TM 




1, 

J 


Q 0347-383 


17.7 


3.220 


800 


1.45 


1.15 


TN 


Jan. 9, 2003 


1 


Q 0420-388 


16.9 


3.120 


600 


1.37 


1.06 


CL, WI 


Dec. 8, 2002 


1 


Q 0913+0715 


17.8 


2.787 


800 


1.45 


1.18 


TN 


Jan. 9, 2003 


2 


Q 1151+0651 


18.1 


2.758 


900 


0.94 


1.17 


TN, CL 


Jan. 25, 2003 


2 


Q 1209+0919 


18.5 


3.291 


1500 


0.78 


1.80 


CL, TN, TK 


Jan. 1, 2003 


2 


Q 1223+1753 


18.1 


2.945 


900 


0.66 


1.35 


TN, CL 


Jan. 25, 2003 


2 


Q 2139-4434 


17.7 


3.214 


800 


0.85 


1.33 


TN 


Apr. 30, 2003 


3 



" Legend: PH-Photometric, CL-Clear, TN-Thin cirrus, TK-Thick cirrus, WI- Windy. 



1: IWorseck et ai] j2007b : PH conditions. 

2: SDSS. 

3: IVeron-Cettv & VeronI i200g) . 
4: Slit loss corrected only. 



mated local continuum and the resulting transmission spectrum 
(see online material for the complete set of quasar spectra). 

3. Analysis 

3.1. The flux transmission technique 

The distribution of absorption lines along a line of sight (LOS) 
towards a quasar is usually expressed as a function of red- 
shift z, column density A^hh ™d Doppler parameter b in the 
form d-'n/ (dzdNwidb) - T](z,Niii,b). Due to our limited spec- 
tral resolution we could not perform single a bsorb e r an alysis. 
We followed instead the approach proposed bv IZuol d 19930 and 
[Liske & Williger (2001) to link the line number density to the 
evolution of the effective optical depth. The resulting evolution 
depends on redshift as 



Teff =B(1+Z) 



y+1 



(1) 



where B and y are sensitive to the resolution and the detectable 
column density range and the observable effective optical depth 
defined as the optical depth at the average transmission over a 
predefined wavelength interval: e"^="^ =< e >. 

In order to account for local fluctuations of the ionising ra- 
diation field, we follow the approach by BD088 which assumes 
intervening absorbers to be in photoionisation equilibrium with 
the local ionising field; furthermore an empty space, and no flux 
attenuation except geometric dilution. The modification intro- 
duced in the optical depth then becomes 



Teff = B(1+Z)''+'(1+W) 



(2) 



where oj is the ratio between the photoionisation rates of the 
quasar and the background and /3 the slope in the column density 
distribution. Assuming a constant UVB in the range of redshifts 
of our sample, an equal spectral energy distribution of QSOs 
and background at v > vll, and pure hydrogen absorbers that 



are isothermal, homogeneous, and randomly distributed along 
the LOS, we find 



U) = 



/v( iLL(l+Zc)) 1 ( dL(z^,0) 

4Kjy (1 +Zc)\dL(Zn,Zc) 



(3) 



with Zc being the cloud redshift, diizq, 0) the luminosity distance 
of the QSO as seen fro m the Earth, and diizg , Zc) as seen from the 
cloud. As clarified bv lLiske & WiUiged (|2001|) . /v(/Ill(1 + Zc)) 
is the flux at the Lyman limit which has to be weighted by a 
bandwidth correction. We computed then the normalised opti- 
cal depth (also called ^) which is the deviation of the detected 
optical depth from the one expected in the Ly Forest 



Teff 



B(l +z)>'+i 



(1 + w) 



(4) 



In order to quantify the reduction of ^ close to a QSO we 
need to know the parameters B and y quantifying the redshift 
evolution of the Lyman forest. These valu es are typi c ally deter- 
mined from high-resolution spectroscopy. Kim et al. ' (2002) ob- 
tained B - 0.0032 and y = 2.37, which we used as starting val- 
ues to compute the normalised optical depths for all QSOs over 
the full spectral range. This resulted in slightly too high average 
^ outside the proximity effect zones, where Eq.[T]should hold and 
produce a mean ^ of unity. We corrected this slight mismatch by 
adjusting B until we reached ^a>^o ~ 1 for the Lyff forest region 
in all spectra, excluding the proximity effect zone. This resulted 
in final normalisation parameter of B = 0.0041 . We will quantify 
the impact of different normalisation strategies when presenting 
the results for the U V background in Sect. 14.11 Finally the slope 
of the column density distribution was set thro ughout the paper 
to /? = 1.5 (e.g. lHu et al.lll995l;lKim et"aLll2002l) . if not expHcitly 
written. 

In order to reveal the proximity effect, we now searched for 
a systematic departure of the normalised optical depths ^ from 
unity for large values of lj. The result for our combined sample 
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Table 2. Redshift estimates with associated errors for the quasar sample resulting from different emission lines ' 



QSO ZLyff 4iii+Oi ^Cii zsiiv+Oiv] zciv ZLit Reference 

or- = 0.003 (T, = 0.003 cr, = 0.005 a-, = 0.003 a- = 0.003 zsiii+Oi - zsiiv+Oiv] 



CTQ 0247 


3.008 


3.025 


3.021 


3.016 


3.008 


0.009 


3.020 


1 


CTQ 1005 


3.201 


3.205 


3.210 


3.196 




0.009 


3.210 


1 


CTQ 0460 


3.134 


3.139 


3.135 


3.128 




0.011 


3.130 


1 


H 0055-2659 


3.650 


3.665 


3.659 








3.655 


2 


HE 0940-1050 


3.081 


3.086 


3.088 


3.074 


3.059 




3.068 


3 


HE 2243-6031 


3.004 


3.010 


3.008 


3.005 


3.004 


0.005 


3.010 


4 


HE 2347-4342 


2.877 


2.885 


2.886 


2.870 


2.861 


0.015 


2.885 


5 


PKS 2126-15 


3.279 


3.285 


3.286 


3.270 




0.015 


3.267 


3 


Q 0000-26 


4.100 












4.098 


6 


Q 0002-422 


2.765 


2.767 


2.765 


2.757 


2.756 


0.010 


2.767 


3 


Q 0347-383 


3.213 


3.220'' 




3.209 






3.222 


12 


Q 0420-388 


3.117 


3.120 


3.120 


3.109 




0.011 


3.110 


2 


Q 0913+0715 


2.784 


2.787 


2.786 


2.779 


2.767 


0.008 


2.785 


7 


Q 1151+0651 


2.755 


2.758 




2.752 


2.754 


0.006 


2.762 


8 


Q 1209+0919 


3.292 


3.291 


3.289 


3.278 




0.013 


3.291 


9 


Q 1223+1753 


2.935 


2.945 


2.944 


2.936 


2.930 


0.009 


2.936 


10 


Q 2139-4434 


3.211 


3.214 


3.210 


3.197 




0.017 


3.230 


11 



t: Wavelength used to esti mate th e redshifts are: Lya = 1215.67 A, Sill + Ol = 1305.77A, Cll = 1335. 30A, Siiv + Oiv] = 1396.76A, 

Civ = 1 549. 06A( Morton 2003h 
a: Taken as systemic redshift. 

b: These redshifts were computed from the average shift between the redshift measurements of the Si n+0 1 and Si iv+0 iv] emission lines. The 

average redshift shi ft is ab out 0. 011. 

Ref : (1) Lopez et all ll2q0lh . (2) lOsmer et alj (Il994h (3^ iRollinde et al.l (l2005h (4) [Lopez et alj (|2002|). (5) iReimers et al.l (Il997h. (6) 
[Schneider et al. (199 ih. (7)|Pettini et al.1 j l997l). (S " ) IVeron-Cettv & Verorj j200a) . (9) IStorrie-Lombardi & Wolfd j2000h . (10) iHewett et al.l 
( Il995h . (11) Hawkins &^rorJ ( 119931) , (12) ISteidell fl99oF ! 



is shown in Fig. |4l while Fig. |5] displays the run of ^ versus co 
for each individual QSO line of sight. Before we discuss these 
results we want to briefly describe our approach to quantify the 
statistical and systematic errors. 

3.2. Error estimates from synthetic spectra 
3.2.1. Method 

Realistic models of the Ly forest, as already developed by 
[Zhang et al. (1995), invoke the baryonic component in CDM 
simulations to trace the absorption line properties along sight 
lines towards QSOs. However, in a first approximation, such dis- 
tributions can be considered as random processes, mathemati- 
cally governed by Poisson statistics. Following this assumption 
we performed extensive Monte-Carlo simulations to study the 
error budget, in particular systematic errors arising from limited 
spectral resolution. 

Each given simulated hne of sight was populated with lines 
distributed using dn /dz oc (1 +zy as line number density distribu- 
tion leading to = B(l +zy'^^. The algorithm continues to add 
absorption features until the effective optical depth r eaches the 
expected value of Tgff using the best fit constraint by Kim et alj 
(2002) B = 0.0032 , y = 2.37. The column density distribution is 
given by /(Nhi) ^ where the slope is ~ 1.5. The Doppler 
parameter distribution is given by dn/db cc b^^exp^-b^-Zb*^ 
where bo- — 24 km/s dKim et alj|2001i) is a parameter depend- 
ing on the average am plitude of the fluctuatio ns in the velocity 
space of the absorbers dHui & Rutledgel[l999l) . Each absorption 
feature was modelled as a Voigt profile and once a transmission 
spectrum was computed, we multiplied it by a QSO template 
spectrum as described in the next section. 




3000 



4000 5000 
Wavelength [A] 



6000 



Fig. 2. Top panel: Average ratio between the fitted and input 
continuum for the five sets of simulated QSOs. Bottom panel: 
Standard deviation profiles relative to the above systematic bias. 
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Fig. 4. Normalised Optical depth versus pro- 
file for the combined sample of 17 quasars, 
binned in steps of A log ai = I. The signature of 
the proximity effect is clearly visible for both 
continuum corrected and uncorrected profiles. 
The curved lines are the best fit of the simple 
photoionisation model to the data, correspond- 
ing to a UV background of log(y,,) = -21.03 
(solid line), in units of [erg cm^-s"' Hz"' sr"'], 
while the green long-dashed line correspond to 
the best fit to the continuum uncorrected data. 
The horizontal line refers to the case of no 
proximity effect. 




Fig. 3. Statistical (shot noise) uncertainties of the normalised op- 
tical depth close to a quasar, estimated from the Monte-Carlo 
simulations. 



3.2.2. Quasar SEDs 

We generated a set of 200 artificial quasar spectral energy distri- 
butions (SEDs) via the principal component method as described 
by Suzuki (2006). Each rest-frame quasar spectrum can be de- 
composed as Xi(/l) - fJ.iA) + S(ci ■ pi{A)) with a mean spectrum 
fi, the principal component spectra p-, and the coefficients Ci. We 
adopted the principal components bv lSuzuki et al.l (l2005h . who 
determined fi and p-, at 1020 A < ^rf < 1600 A from 50 HST 
FOS spectra of low-redshift quasa rs. The coeffi cients Ci are ap- 
proximately Gaussian distributed (Suzuki 2006). After generat- 
ing the 200 SEDs by drawing the Ci from their Gaussian distri- 
butions, we convolved it with the instrumental profile and added 
random Gaussian noise to reproduce our observations. 

In order to investigate how different QSO emission redshifts 
affected the error budget, we simulated quasars at five typical 
redshifts z = 2.0, 2.5, 3.0, 3.5, 4.0, (denoted in Figs. HHD with 
numbers 1, 2, 3, 4 and 5, respectively) and luminosities repre- 
sentative for our sample. 

We then normalised our spectra as performed for the actual 
observations. These simulations created a database to quantita- 
tively assess the two dominant sources of error: the limited num- 
ber of absorbers per individual line of sight (cosmic variance), 
and the misplacement of the continuum because of line crowd- 
ing. In the following we consider both eiTor sources in turns. 



3.2.3. Systematic errors: Continuum placement 

For each simulated spectrum we computed the ratio between the 
fitted and the input continuum. We then averaged the ratio over 
all realisations at each emission redshift. The result is shown in 
the top panel of Fig.|2l for the five redshifts adopted. Expectedly, 
the highest deviation is always in the wings of the Lya emission 
line. It can also be seen that the increasing line number density 
with redshift causes a gradually growing systematic error We 
use this profiles to correct our automatic continuum estimates 
and use the standard deviations shown in the bottom panel of 
Fig.|2]as contribution to the uncertainties. The effect of this cor- 
rection can be seen in all the proximity effect plots as difference 
between the green triangles and the black dots in Fig.|4]|5] 

We note in passing that a second source of systematic errors 
would be the presence of metal line systems in the proximity 
effect zone. Our spectral resolution is insufficient to identify in- 
dividual metal lines in the Lyman forest, and any such absorp- 
tion present, but unaccounted for, will systematically increase 
the normalised optical depth ^ and thus tend to mask the prox- 
imity effect. 

3.2.4. Statistical errors 

We modelled the statistical error of the measured optical depth 
along individual lines of sight as arising from Poissonian shot 
noise due to the limited number of absorbers in each simulated 
spectrum. For each stack of simulations we computed the mean 
and standard deviation of ^ per log a> bin. The results are shown 
in Fig. [3] which demonstrates that the standard deviations cr^ 
are always considerably bigger than the continuum dispersion. 
However, recall that we fully account for both random and sys- 
tematic errors. As expected, the statistical errors become bigger 
towards lower redshifts due to the smaller line number density. 
The simulation results were then used to describe cr^icj, z) with 
a simple polynomial parameterisation. Without introducing the 
proximity effect, the statistical error at high a> would be much 
larger due to cosmic variance on very small scale. 

In order to estimate the statistical eiTor for the combined 
analysis of the full sample we ran a new set of simulations. Here 
we generated 10 random lines of sight with the emission redshift 
of each of the 17 quasars in the sample, and computed the sta- 
tistical scatter after averaging over the 17 contributions to each 
cj value. We found cr^ to be essentially independent of oj in this 
case. 
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4. The proximity effect in the combined sample 

4.1. A new estimate of the UV background 

Figure |4] summarises the main results regarding the combined 
sample of 17 QSO spectra. We see a highly significant reduction 
of normalised optical depths at log w <; 0, i.e. for the zone where 
photoionisation due to the local quasar-induced radiation field 
is expected to prevail over the metagalactic UV background. As 
already demonstrated by BD088, this turnover can be used to 
constrain the mean intensity of the UVB. We adopted the fitting 
formula 

F(co) = (l + ^) (5) 

with a being the free parameter and w_2i is the value of oj rela- 
tive to a reference value of the UVB, J*^^ = 10"^' erg cm^^s"' 
Hz ' sr '. We then applied a straightforward minimisation to 
search for the best-fitting value of a = J(yi±)l J*2y 

In doing this computation we found that the bin size A log w 
has some moderate effect on the resulting best-fit value of the 
UVB. If the data are merged with very small or even without 
any binning, a will be biased towards low values because of the 
substantial scatter in ^ at very small w (for log(w) % -1.5), due 
to the strong effects of shot noise for this narrow log(w) interval. 
On the other hand, too large bins will tend to hide the signature 
of the proximity effect, thus make the UVB appear stronger than 
it really is. As a reasonable compromise we chose Alogw - 
1, upon which also Fig. |4]is based. An additional effect which 
tend to change the estimation of the UVB is the normalisation 
used to compute ^ (see eq. |4|i. We address this pro blem with 
two di fferent strategies. We use the normalisation by iKim et al.l 
(120021) and our normaHsation (B - 0.0041) to reach ^^j^q ~ 1 
for the combined set of Lya forest regions. Tab. |3] summarises 
our results revealing a maximal dispersion of about ~ 0.9. We 
decide to use B - 0.0041 for our results. 

As best fit value we obtain J{vi±) = 9 ± 4 x 10"^^ erg cm"^ 
s"' Hz"' sr"', or in logarithmic units, logy(vLL) = -21.03^q22. 
Using a slightly narrower bin size of A log u - 0.7 lowers /(vll) 
by about 0.05 dex. 

This estimate of the UVB intensity is in very good agree- 
ment with all recent measurements base d on a wide ra nge of 
techniques and data sets. For example, Scott et al.l (12000) ob- 
tained /(vll)/-/*2i = 0.7^Q44, applying line count statistics on 
more than hundred spectra at ~ lA resolution. More similar to 
our approach. Liske & Williger (2001) used the flux transmis- 
sion statistic on 10 QSO spectra with ~ 2 A resolution and a 
S/N of ~ 40, obtaining 7(vll)//*2i = 0-35!o 13- Not much has 
yet be en published using high resolution spectra. Giallongo et al. 
(il996) obtained y(vLL)//*2i = 0.5 + 0.1 and lCooke et al.l (119971) 
J{vi^i} / J*2i - 0.8^y 4, again close to our value even though they 
are all smaller (up to a factor of about 2). This order of magni- 
tude is also con sistent with predictions based on the QSO lumi- 
nosity function (iHaardt & Madau,1996.) . 

4.2. Dependence on model parameters 

Concerning the evolution of Tetf in the Ly forest (see eq. [TJ, we 
regard only the normalisation B as variable. The slope 7 has 
been estimated by s everal authors with much higher accuracy i n 
the past years (e.g.. iKim et ani2001l llOOl ISchave et al.ll2003h . 
We consider only two norma lisations: a conservative one with 
B = 0.0032 dKim et al.ll2002l) . which returns a higher UVB es- 
timate since ^^^0 ^ 1, and B = 0.0041 which we adopt since it 



Table 3. Estimations of the UV background /(vll)/-/*2i (com- 
puted in units of erg cm"^s"' Hz"') from different model param- 
eters (fi and B) and binsizes A log oj. 





B = 


0.0032 


B = 


0.0041 


Alog(w) 


0.7 


1.0 


0.7 


1.0 


yS= 1.4 


1.01 


1.12 


0.40 


0.51 


/?= 1.5 


1.89 


2.04 


0.88 


0.93 


13= 1.6 


2.83 


3.01 


1.19 


1.41 



leads to ^^^o - 1- We believe that the discrepa ncy between our 
adopted value and the value bv lKim et al. I (l2002l) is due to a com- 
bination of line blending and resolution effects. Our results are 
sensitive to the slope of the column density distribution /? which 
sets the steepness of ^{oj) in the BD088 ionisation model. The 
column density distribution is very well a pproximated by a sin - 
gle power law with yS - 1.5 over 10 dex dPetitiean et al.|[T993l) . 
However, it has been shown that the slope changes somewhat 
with redshift (Kimet al. 2002). Since the lines are unresolved 
in our low-resolution spectra, we fix a single power law distri- 
bution and vary its assumed value by Ay6 = 0.1 to estimate the 
robustness of our UVB measurements. The Doppler parameter 
distribution has no direct impact on our results since the BD088 
model simply assumes an isothermal H i distribution. 
Table |3] summarises the dependence of our UVB estimates on 
model parameters. There is a substantial scatter in the estima- 
tions and both B and /3 play a central role. The binsize has only 
a small (but still detectable) effect. 

5. The proximity effect in individual lines of sight 

5.1. Results 

The proximity effect is generally seen as a statistical phe- 
nomenon, which may or may not be detectable in individual 
spectra. The high S/N of our spectra motivated us to search for 
proximity effect signatures in each of our 17 quasar spectra. The 
basic approach was essentially identical to that of the combined 
analysis. Compute the mean normalised optical depth for a given 
line of sight within a given bin A log co and check whether ^ sys- 
tematically decreases for large values of oj. We set the same nor- 
malisation for every object as for the combined analysis. The 
o) scale was now fixed by assuming the value of the mean UV 
background intensity from the combined analysis. The results 
are displayed in Fig. |5] one small panel per quasar The error 
bars are now of course dominated by Poissonian shot noise, es- 
timated from the simulations as described above. In each panel 
the expected run of ^ with oj is shown as the dotted red line, as- 
suming that the metagalactic UV background is constant and has 
the same spectral shape as the quasar. We also show the profile 
before the systematic continuum correction and its best fit model 
(green triangles and line). 

Figure|5]demonstrates that in all except one case, ^ decreases 
substantially from left to right. Thus we can say that the prox- 
imity effect is detected in 16 out of our 17 quasar spectra. In the 
majority of spectra, the ^-cj profile is even formally consistent 
with the prediction based on the combined analysis. In a number 
of cases the data seem to be (mildly) discrepant with the predic- 
tion; we discuss some of these cases below. 

We also applied the above fitting procedure to each spec- 
trum separately; the results of that exercise are shown as the solid 
curves in Fig.|5] Table |4] summarises the fit results. Apparently 
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Fig. 5. Search for proximity effect signatures in individual lines of sight, for all 17 QSO spectra. Each panel shows the normalised 
optical depth ^ versus oj in the same way as Fig. ID with the best-fit model of the combined analysis superimposed as dotted red 
lines. The solid lines delineate the best fit to each individual QSO as described in the text. For CTQ 0247 we also plot the best fit 
excluding the strong absorption at z ~ 3.014 (purple curve). The panels are sorted in order of decreasing strength of the proximity 
effect (horizontal displacement of solid and dotted lines). HE 2347-4342 has no detectable proximity effect in H i. 
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the value of the fitting parameter a shows significant scatter be- 
tween the different quasars. The value of a describes the horizon- 
tal offset of the solid curve relative to the dotted curve and in the 
following is regarded as a quantitative measure of the strength 
of the proximity effect signal (in the sense that a large a means 
a weak proximity effect). This does of course not imply that the 
UV background fluctuates by a similar amount. While statistical 
errors certainly contribute to the scatter, one may also reinter- 
pret the parameter a as a measure of the flux of the quasar at the 
Lyman limit. This flux may even not always have been constant 
over the light travel time across the proximity effect zone; we 
return to that point in Sect. l5.4l below. 

We now briefly comment on two objects where the proxim- 
ity effect appears to be extremely weak or absent. In the case 
of CTQ 0247, there is a strong associated absorption system at 
Zabs - 3.014, corresponding to logw < 1. Removing this ab- 
sorption manually from the spectrum and redoing the analysis 
yielded a proximity eflFect signal perfectly consistent with the 
prediction from the combined analysis. 

5.2. A hidden proximity effect for HE 2347-4342? 

For HE 2347-4342, we see no evidence at all of a downturn of 
^(oj). The absence of any proximity effect was noticed already by 
iReimers et al.l (Il997h upon mere visual inspection of the spec- 
trum. Again there is a conspicuous strong associated absorption 
system. Although in this case a removal of that system does not 
dramatically improve the proximity effect signal, the absorber 
may nevertheless play an important role in this line of sight, as 
it may attenuate the ionising radiation field towards lower red- 
shifts. This is supported by the following simple calculation. 

The total measured H i column den sity in the associated sys- 
tem is of the order of 2 x 10 cm" jFechner et al.ll200l . We 
assume this to be located in an absorbing slab of matter very 
close to the QSO. Knowing the Lyman limit flux of the QSO, we 
can predict the value of oj at the location of the slab, log w ^ 2. 1 ; 
this immediately relates to a predicted reduction of H i column 
density in the slab of ~ 2 dex. Thus, the column density of the 
same absorber without the QSO ionising radiation would be of 
the order of 10'^ cm"^. This would be sufficient to render a re- 
maining proximity effect for the line of sight undetectable. 

We note that Fechner et alJ ("2004) detected evidence of a 
hard radiation field from a detailed photoionisation modelling 
analysis of metal absorption lines in this system. A similar trend 
is apparent in our recent investigation of the He ii Lyman forest, 
comb ining VLT and F USE high-resolution spectra of this quasar 
IWorseck et al.l (l2007h . Thus, while the traditional Hi proximity 
effect is clearly absent in HE 2347-4342, most probably due to 
excess absorption close to the QSO, ther e are clear signs of a 
'proximity effect in spectral hardness' (cf. IWorseck & Wisotzkil 
[2OO6) for this object. We therefore conclude that all 17 quasars 
in our sample show evidence of a genuine proximity effect. 

5.3. Variations in tlie strengtli of tlie proximity effect 

We now investigate how much of the scatter in the fitting pa- 
rameter a might be attributed to uncertainties, or shot noise, or 
intrinsic dispersion of other relevant properties. 

Figure |6] shows that formally the best-fit values of a, be- 
fore and after continuum correction (green and black histograms 
respectively), extend over several orders of magnitudes, with a 
standard deviation of criog(^) ~ 1 (black histogram). 



Table 4. Input and results used for the photoionisation model 
fits, /ll is the predicted Lyman limit flux of each quasar. In 
the treatment of individual lines of sight, log(fl) is the fitted pa- 
rameter used to quantify the strength of the proximity effect. 
CTQ 0247 is listed twice, with and without the associated ab- 
sorption system. 
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The large dispersion can to a small part be explained by un- 
certainties in quasar magnitudes and redshifts, which both affect 
the computation of a>. For those objects where we have rather 
accurate V magnitudes (with errors < 0.1 mag), the flux scale 
is accurate to within 9%. For the remaining objects the flux un- 
certainties might be as large as 30%. The direct effect on log(fl) 
is an uncertainties of 0.03-0.1 dex. In addition, redshift errors 
(Tab. [TJ might shift log(a) by 0.05-0. 1 dex. We conclude that 
these uncertainties cannot be the main source of spread in the 
estimated log(a) values. 

We employed our Monte-Carlo simulations to estimate the 
expected scatter solely due to statistical shot noise errors (i.e., 
cosmic variance). To this effect we first systematically reduced 
the optical depths in the simulated spectra following strictly the 
theoretical proximity effect prescription (Eq. |4]i, after which we 
'remeasured' the proximity effect and its strength a by fitting the 
artificial data in the same way as the observed ones. The result- 
ing histogram of a values is superimposed in Fig.|6j the distri- 
bution is approximately Gaussian with a standard deviation of 
o'iog(fl) = 0.64. This is substantial and implies that 14 our of 17 
of our quasars are located within +2cr expected for pure random 
errors. Nevertheless, 3 quasars from our sample are located out- 
side of this interval, which is a bit much to declare them all as 
outliers. 

This means that other effects play a role, and that these ef- 
fects may lead to gross deviations from the simple expectation. 
Possible mechanisms might be, for example, a dispersion in 
spectral indices for the quasars, long time scale variability of 
the quasars, or a strongly fluctuating UV background. We first 
consider the effects of non-uniform spectral indices. 

In the standard analysis of the proximity effect as introduced 
by BD088, the assumption of a uniform quasar spectral index 
is obviously wrong, but the averaging over samples of quasars 
makes the analysis very insensitive to any intrinsic dispersion. 
When considering individual quasar lines of sight as we do here, 
this assumption may be more harmful. We investigated this by 



A. Dall'Aglio et al.: The line-of-sight proximity effect in individual quasar spectra 



9 




Fig. 6. Distribution of the parameter log(fl) for the fits to the 17 
individual spectra with and without continuum correction (black 
and green hashed histogram respectively). Overplotted is the dis- 
tribution resulting only from shot noise in the simulated spectra 
(solid line), and the expected distribution of log(a) due to disper- 
sion of spectral indices (blue gaussian). 



extending the original calculation by BD088, solving the exact 
integrals for the ratio of the photoionisation rates of the quasar 
and the background and using the photoionisation cross section 
for a pure hydrogen cloud, so that 



hv 



r 



4;ry°(v) 
hv 



o-(v)dv ^ 



cr (v) dv 



a — 3 



(6) 



where Woid is the know expression of eq. [3] and k ~ -1.8 
is the b ackground spectral index following the model by 
iHaardt & M adau C1996) . Expressing everything in logarithmic 
units we arrive at 



log(w) = log(Wold) + Criog(a) 



(7) 



Real quasars show a dispersion of spectral slopes of era = 0.25 
around ff = -0.46 (derived from the SDSS Q SO composite spec- 
trum and its standard deviation published by I Vanden Berk et al.l 
1200 li) . The values of oj are therefore offset (i) by Alog(aj) - 
0.142 on average, and (ii) by a random component with a stan- 
dard deviation of criog(a) = 0.06 (included as the narrow Gaussian 
in Fig. |6]l. Obviously, the dispersion in spectral slopes is too 
small, by a long way, to account for the observed distribution 
of log(fl). 

Among the two other mentioned options, intrinsic fluctua- 
tions of the UV background might also contribute, but numeri- 
cal simulations suggest that substantial fluctuations are only ex- 
pected for redshifts higher th an those covered here jCroft et al.l 
[l999tlMcDonald et al.ll2005l) . 

The presence of metal transitions falling in the proximity ef- 
fect zone or a still imperfect continuum placement might also 
enhance the scatter. 

A very plausible eff'ect, on the other hand, would be signif- 
icant QSO variability over timescales of the light travel time 
across the proximity eff'ect zones which is of the order of 10^ 
years; there is no reason to expect that quasars always maintain 
their radiative output over such long periods. If they do not, then 
the observed Lyman limit luminosity will not be the same as that 
received by clouds along various points along the line of sight. 
We will discuss a specific aspect of this effect in the next section. 
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Fig. 7. Top panel: Expected ^(oj) model for different finite life- 
times of the quasar. The solid curve represents the BD088 model 
of the proximity effect, with a QSO of infinite lifetime. Bottom 
panel: Tentative application of the above model to the QSO 
Q 0002-422. The dotted line is the predicted ^{a>) curve for 
infinite lifetime, the dot-dashed line delineates the curve for a 
fiducial finite lifetime of 68 Myr. 



5.4. Finite quasar age 

The simple theory of the proximity effect implicitly includes the 
assumption that quasars shine for an infinite time, or at any rate 
for much longer than the light crossing time of the proximity ef- 
fect zone. This may be wrong, and we now ask specifically what 
would happen if a quasar is abruptly switched on within, say, less 
than a few Myrs before the observation. As soon as the quasar 
starts to radiate in the UV, an over-ionised "sphere" will start to 
expand around the quasar (always assuming spherical symme- 
try). The new equilibrium state with ionising photons from the 
local source and the UV background - requires at least some 
~ lO'* yrs to establish, which means that for even younger QSOs 
the proximity eff'ect will be absent. 

The light travel time between two redshifts along a single 
line of sight can be calculated for a A-Universe as 



t(z) 



3Ho va 



:asinh • 



(8) 



(iPeacockl 11999 ^. allowing us to convert redshift intervals into 
light travel time diff'erences. Switching a quasar on can be ex- 
pressed in the model by introducing a step function into the oj 
profile. 



r 

(^(Z) = { /v(iLL(l+-- 

4;ry, 



+--c)) 1 / fl'z.fc|,0) f 



for Z < Ziife 
forziife < z < 



(9) 



Far away from the QSO (low cj values), the QSO has no effect 
yet. When passing into the over-ionised region, the theoretical 
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profile ^(lo) drops abruptly, joining the BD088 prescription of 
the proximity effect. This is shown in the top panel of Fig. [T] 
where we present the expected profiles for different assumed 
ages. For presentation purposes, this model was computed for a 
QSO at zq = 2.76 and /v(/Ill(1 + Zq)) = 3.46 ■ lO^^? erg cm-^g-' 
Hz 1. 

It turns out that observable signatures can be expected for an 
interesting range of ages. Up to fQso ~ 4 Myr the turnover hap- 
pens so close to the quasar that even for a very luminous QSO it 
will be hard to detect, given the inevitable shot noise limitations. 
Above of at most ~ 500 Myr, on the other hand, the differences 
between the models with and without lifetime will wash out 
completely because ^(co) is expected to be close to unity anyway. 
The range between these extremes, 4Myrs ^ fgso ^ 400Myrs, is 
interestingly close to quasar lifetimes estimated by other, usually 
muc h more in direct methods, such as modelling of QSO accre- 
tion dYu & Tr emaine 2002; Hopkins et al. 2006) or analysis of 
QSO clustering (Croom et al. 2005). It might thus be possible to 
detect age effects of quasars by studying their proximity effect 
signature. 

We searched our quasars for possible examples of such a 
step feature, and found one (highly tentative) possible exam- 
ple: Q 0002-422, shown in the bottom panel of Fig. |2l seems 
to match the expected pattern for a 'recently born quasar' quite 
well. The normalised optical depth ^ remains constant at ~ unity 
over the Lyman forest until log(w) - -1.5 is reached, where it 
rather abruptly drops to much lower values. We chose a smaller 
binning in log co for this plot, compared to the previous figures, 
in order to highlight the relatively sudden drop. Applying our 
simple model we obtain an age of Iqsq ^ 68;^^' Myrs, where 
the errors are based on simply assuming the uncertainty to be 
+0.35 bin in log a>. Of course we do not claim to have measured 
the age of that quasar with 20 % accuracy. The difference in the 
goodness-of-fit between finite and infinite quasar ages is not for- 
mally significant. But the exercise shows that it is possible to 
derive empirical constraints on quasar ages from proximity ef- 
fect signatures, which might even be enhanced by using data of 
higher spectral resolution. 



6. Conclusions 

We have presented new evidence of the line-of-sight proxim- 
ity effect as a universal phenomenon occurring in the spectra 
of high-redshift quasars. Even though our spectra are limited 
in spectral resolution, their high S/N and the power of the flux 
transmission method has enabled us to demonstrate the presence 
of the effect for every single line of sight, for the first time. 

Our estimate of the mean UV background intensity for the 
redshift range 2.7 < z < 4 is log(7v) = -21 .03;^q 22, in very good 
agreement with literature values for similar redshift ranges. We 
made a careful assessment of the error budget using extensive 
Monte-Carlo simulations. The errors are clearly dominated by 
cosmic variance, which implies that better spectral resolution 
would not necessarily have a dramatic impact on the measure- 
ment accuracy. Of course, high resolution spectra would be ad- 
vantageous for a more detailed analysis of several other aspects 
of the proximity effect, such as the effects of gravitation al clus- 
tering of absorbers near the QSO (iRoUinde et al.ll2005h . which 
may overe stimate the value of the U V background by up to a 
factor of 3 dLoeb & Eisensteinlll995h . 

We have quantified the strength of the proximity effect in in- 
dividual spectra and find that this shows a higher dispersion than 
expected from only statistical shot noise errors. Among the most 



likely contributors for this additional dispersion are again gravi- 
tational clustering of absorbers near the QSO, or QSO variabil- 
ity over very long timescales; fluctuations of the UV background 
are also possible but unlikely to play a major role at this redshift 
range. We presented a speculative, but conceptually simple ob- 
servational test to search for signatures of finite quasar ages in 
the optical depth profiles derived for a single quasar, and we even 
tentatively identified a candidate where such a pattern might be 
visible. 
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